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Abstract 

We investigate through analytic calculations the surface plasmon dispersion relation for mono- 



layer graphene sheets and a separated parallel pair of graphene monolayers. An approximate 
form for the dispersion relation for the monolayer case was derived, which was shown to be 
highly accurate and offers intuition to the properties of the supported plasmon mode. For parallel 

graphene pairs separated by small gaps, the dispersion relation of the surface plasmon splits into 

D 

two branches, one with a symmetric and the other with an antisymmetric magnetic field across 

the gap. For the symmetric (magnetic field) branch, the confinement may be improved at re- 
duced absorption loss over a wide spectrum, unlike conventional SP modes supported on metallic 

surfaces that are subjected to the trade-off between loss and confinement. This symmetric mode 

O ■ 

Q ■ becomes strongly suppressed for very small separations however. On the other hand, its antisym- 



metric counterpart exhibits reduced absorption loss for very small separations or long wavelengths, 
serving as a complement to the symmetric branch. Our results suggest that graphene plasmon 



> 

OO 

£T) ' structures could be promising for waveguiding and sensing applications in the mid-infrared and 



terahertz frequencies. 

PACS numbers: 73.20.Mf,73.25.+i,78.66.-w 



I. INTRODUCTION 

Surface plasmonsi (SPs) are charge density waves propagating at the interface between a 
conductor and a dielectric medium. SPs are transverse-magnetic (TM) polarized with modal 
fields that peak at the interface and decay exponentially away into both media. Noble met- 
als such as gold and silver have traditionally been the material of choice to support SPs, 
providing for lossy but well-confined plasmon fields at visible frequencies. To achieve SP res- 
onance in the near-infrared regime and beyond, metallic films corrugated with nanoholes^ or 
nanorods,- and multi-layered metmaterial structures^ have been proposed. More recently 
graphene, a single layer of carbon atoms gathered in honeycomb lattice, has also been con- 
sidered for supporting surface plasmon modes^— at infrared and terahertz frequencies. The 
transport properties of graphene can be readily tuned by application of a gate voltage^ 1 ^, 
opening the possibilities to various plasmonic devices^2r— . Moreover, several techniques 
have been proposed to experimentally realize the excitation of plasmons on graphene^ 1 ^. 

The SP dispersion of isolated monolayer graphene (~ 0.34 nm thick) and monolayer 
graphene deposited on dielectric substrates have been studied by several authors.— ~— The 
plasmon dispersion for parallel graphene pairs,— >^ multilayer graphene stack,— and interca- 
lated graphite^ have also been theoretically investigated. For a parallel pair, the geometry 
approaches the bilayer graphene when their separation is small (few angstroms), and effects 
of interlayer hopping could become important^ 1 ^. 

The transport properties of graphene can be modeled with the linear response the- 
ory(Kubo formula) I&2&2L22, In the nonlocal random phase approximation (RPA), analytical 
expressions that are valid for arbitrary frequency u and plasmon propagation wave vector 
j3 may be derived in the limit of zero temperature £i22i2L2iL2£ Alternatively in the local limit 
where spatial dispersion effects are negligible,— i.e., u » r~ l and u > (3vp, analytical 
expressions for the optical conductivity (a) of graphene can be obtained from the Kubo 
formula for finite temperatures T. Here r is a phenomenological electron relaxation time, 
vf ~ 10 6 m/s is the Fermi velocity, and — |/3|. As dynamic control of the SP disper- 
sion holds potential promise for the design of various active devices such as optical switches 
and modulators, it is important to study the characteristics of the SP modes such as their 
confinement and spatial distribution of dominant field components, so as to gain further un- 
derstanding of their physical properties. In this work, further to analyzing the SP dispersion, 



we describe within the framework of electromagnetic waveguide theory, how the properties of 
the SP modes supported on monolayer graphene and parallel graphene sheets may be tuned 
over a broad range of frequencies for potential waveguide applications. The characteristics 
of the SP modes are analyzed through their associated electromagnetic mode profiles. In 
addition, we introduce a figure-of- merit (FOM), which in our opinion represents an optimal 
measure that takes into account the compromise between confinement and propagation loss 
of the supported SP modes to quantify the performance of the waveguide. 

To describe the transport properties of graphene in our study, we calculate the local 
conductivity of graphene from the Kubo formula. Let us define the effective index of the 
supported plasmon mode as n = /3/ko, where k^ = 27r/A and A is the free-space wavelength. 
In the local limit (a; >> r^" 1 ), the condition fivp < u translates to n < c/v B w 300 (c being 
the speed of light in vacuum), and the expression for a is2&2& 
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where F(e) = (1 + exp[(e — fi^/KsT])" 1 is the Fermi-Dirac distribution with [i c as the 
chemical potential. In Eq. (JTJ, the first term corresponds to contributions from intraband 
electron-photon scattering and the second term arises from contributions due to direct in- 
terband electron transitions. Integration of the first term leads to^ 

^J^f^^cosh^)], (2a) 

(p c »K B T), (2b) 
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For the interband term, we obtain, after evaluating the numerator of the integrand 
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where G(£) = sinh^/[cosh(/x c /A'sT) + cosh^], and ^ = e/K B T. For finite values of r, 
the integrand in Eq. f l3a|) has no poles along the real £— axis, and is suitable for numerical 



integration. For fi c » KbT, Eq. ( 12b|) shows that cr mtra is directly proportional to the 
chemical potential fi c , while Eq. (13bj) shows that <j mter diverges logarithmically for fuv ~ 2/x c . 
Furthermore, Re(<7 mter ) jumps by an amount of <7o at the interband threshold hu mtcr = 2/i c , 
evident of strong absorption loss for u > cj inter . At high frequencies, u » (fi c ,KsT), the 
graphene conductivity a approaches a constant a rs e 2 /4h, which is close to the measured 
minimum conductivity value of J 



2-nh 



For all following calculations, it is taken that T = 300 K (room temperature, KbT ~ 
26 meV), and that /x c >> KbT . In this limit, Eq. ( 12b|) and Eq. f ]3b|) for the optical con- 
ductivity of graphene holds, and the relationship between the chemical potential \x c and the 
carrier density n c on the monolayer graphene is 

Ai c = \/(hv F ) 2 n c Ti - (irK B T) 2 /3 « hv F ^/n~n . (4) 

Since carrier density of up to 10 14 cm -2 has been realized in experiments^^ 2 -, we set the 
range of chemical potential to 0.2 eV < /x c < 1.2 eV. Based on measured values of the carrier 
mobility^ 1 - 1 ^^, which ranges from 10 3 — 10 4 cm 2 V -1 s" 1 , the relaxation time r is taken to 
be in the order of ~ 10~ 13 s. In graphene, effects of free carrier absorption^ 3 - must be taken 
into account for photon energies above the optical phonon energy 3 ^ hu op = 0.196 eV(A op ~ 
6.3 (Urn) but below hw mtcr . For c<; mtGr > ui > u op , conservation of momentum for the free 
carrier absorption process can be satisfied with the emission of an optical phonon, leading 
to a significant decrease^ in the relaxation time r. A faster relaxation time corresponds to 
higher absorption losses, i.e., an increase in Re(er). The analytical expressions above for 
a (Eqs. ([I])-®) hold for electron and hole bands which exhibit a linear energy dispersion 
relation near the zero bandgap points in graphene^. The range of energy for which this 
linear dispersion relation is valid extends well into the visible spectrum, as shown by a 
recent experiment^. To stay within the regime for which our analysis is valid, we limit 
the calculation of the graphene conductivity a with Eqs. (J2J) and ([3]) to wavelengths in the 
range 0.8 /im < A < 50 /xm (near-infrared to terahertz regime). For long wavelengths, 
the conductivity is dominated by the intraband term, with the real and imaginary part 
scaling as A 2 and A, respectively, in accordance with Eq. fl2bl) . The conductivity decreases 
with wavelength until it approaches the interband threshold, where Re(cr) — y £r , and 
Im(cr) becomes negative. When Im(cr) is negative, TE (transverse electric)-polarized surface 
mode*^£ can be supported on graphene. This mode, whose dispersion lies close to the 



light line, will not be further discussed here. Instead we focus on the TM plasmon modes, 
which offers the potential for smaller effective mode areas compared to the noble metals, 
especially in the long wavelength regime where SPs supported on metals are known to be 
poorly confined. 

The remainder of this paper is organized as follows. In Sec. II, the simplest case where 
monolayer graphene is sandwiched between two dielectric half-space (61,62) is re- visited. It 
is shown that only the antisymmetric mode (labeled as the A mode) is supported in the 
case when e\ and e 2 are lossless dielectrics. By antisymmetric, we mean that the tangential 
magnetic field component exhibits a zero across the monolayer graphene. An approximate 
form for the dispersion relation for the monolayer case is derived, which is shown to be highly 
accurate and offers intuition to the properties of the supported plasmon mode. The case of 
two layers of monolayer graphene with identical chemical potential (/ii = /i 2 ) separated by a 
thin dielectric medium (ed) of thickness d is investigated in Sec. Ill, where the separation d 
is exploited as an additional degree of freedom to tune the properties of the SP modes. For 
the fully symmetric geometry (/i! = /i 2 , ei = e 2 ) and when {e l7 e 2 , e^} are purely real, the SP 
dispersion splits into a symmetric AS branch, and an antisymmetric AA branch, where the 
first A in each label refers to the antisymmetric magnetic field distribution across each mono- 
layer graphene. The second label S (A) refers to the mode with symmetric (antisymmetric) 
magnetic field across the dielectric medium separating the graphene pair. The variation of 
the SP dispersion with wavelength (A), chemical potential (/i c ), and the dielectric constant 
of the separating medium (e^), will be examined. The influence on the dispersion for an 
asymmetric geometry where the chemical potential on each graphene monolayer is different, 
will also be discussed. Finally, we offer concluding remarks in Sec. IV. 

II. SP MODES ON MONOLAYER GRAPHENE 

To begin, we analyze the case where monolayer graphene (at x = 0) is sandwiched between 
two dielectric media (ei,e 2 ) as shown in Fig. 1(a). For a guided wave with wavenumber /3 
propagating along the z axis, the tangential component of if -field along the y axis can be 
expressed as 

, N Hi exp(«ix + \Bz) x < 
H =-- { y ' (5) 

N H 2 exp(— a 2 x + ij3z) x > 



where N is a normalization constant such that the mode has unity Poynting flux at z — 0, 
and 

a 2 = f3 2 -e 3 k 2 , j = {1,2}. (6) 

Here we suppress the time dependence exp(— lout) for notation brevity. The electric field 
components, evaluated from Maxwell equations are 

E Xj = JLh, , (7a) 

U!€j€q 

At x = 0, the continuity of the tangential electric field (E Zl = E Z2 ) yields 

H 2 = -^H ± , (8) 

a 2 ei 

and the boundary condition for the tangential magnetic field, i.e., x x (H 2 — Hi) = erE (x = 

unit vector along x) leads to H 2 — Hi = crE z , or 

H 2 = (1 + ^^-)H X . (9) 

weie 

On substituting Eq. (jSJ) into Eq. (jSJ), we arrive at the dispersion relation 

^ + ^ + ^ = 0. (10) 

The dispersion relation fllOl) bears a striking resemblance with that for the SPs at a metal- 
dielectric interface. If either one of the dielectric media is a metal, then it is clear that the 
SP mode is almost identical to that for the metal-dielectric interface since \a\ lies in the 
range 10~ 4 — 10~ 3 S and u is in the range 10 13 — 10 15 s" 1 . We have verified this fact numer- 
ically, which is consistent with experimental observations where carbon-on-gold substrates 
are used to improve the stability in SP resonance detection of DNA arrays^. Because 
of the similarity to the classical metal-dielectric interface, the metal-graphene-insulator 
case is not interesting in the present context and we shall focus only on the insulator- 
graphene-insulator (I-G-I) structure. To aid in the analysis, we note that the intraband 
conductivity of graphene takes the form similar to the Drude model for free electrons, i.e., 
a intra _ i e 2 Vpy jf[~ c / [y/rfh^u + ir" 1 )], and define an effective plasmon frequency 

r ..(rff)N2 = Vn~cV F e 2 
{ p } ~ e hA 

= -^~T (11) 

e m c A 



where m c = hy/n^f/vp is the effective mass of the Dirac fermions^, and A is the effective 
thickness of the monolayer graphene. For graphene, uip serves as an estimate of the 
plasmon resonance frequency. As opposed to metals with parabolic electron dispersion where 
the plasmon frequency is proportional to nj , Up oc nj for graphene.— >2i Let us note 
that A in Eq. (II ip is not necessarily the physical thickness of the graphene layer, rather it 
is the effective thickness that would allow us to recover the correct dispersion relation when 
the monolayer graphene is treated as a thin metal film or screening layer. To estimate u>p , 
we take A ~ 0.34 nm, which gives \ p e = 2irc/(jjp ~ 300 nm for 0.4 eV < /x c < 1.2 eV. 

Far from the plasmon resonance frequency, the loss of the bound modes or equivalently 
the imaginary part of the wavenumber (fl = 0' + i(3") is small. Consequently, Eq. ([6]) implies 
that the imaginary part of a is small compared to the real part (a = a'+ia"), and a ~ a'. As 
a' must be positive for the bound modes, Eq. (JSJ) implies that the if-field must change sign 
across the monolayer graphene when both the adjacent media are lossless dielectrics, and 
hence only the antisymmetric mode (labeled A) is supported in the case of the I-G-I. Taking 
e 2 = 1 and \x c = 0.8 eV, the real part of the SP mode effective index (n = n' + in" = /3/ko) 
for ei = 1 and 4 is shown in Fig. 1(b). Following the derived results of Ref. [21], we have also 
calculated the SP dispersion relation with the nonlocal RPA, which shows a good agreement 
with Eq. ( TTUj) (see black-dashed curves in Fig. 1(b)). Similar to the antisymmetric SP modes 
supported on thin metal films, it is seen that n' > {y^el, y^ }• As A — > X p e , n! increases 
rapidly. The electromagnetic field components for {ei, e 2 } = {1,1} and {ei, e 2 } = {4, 1}, are 
shown in figs. 1(c) and 1(d) respectively, taking fi c = 0.8 eV and A = 8 fim as an example. 
It is evident that the if -field of the A mode switches sign across the monolayer graphene, 
and has a maximum at the side of the dielectric with a higher permittivity (if e\ ^ e 2 ). It is 
worth noting that because /3 is relatively large compared to ko, the magnitude of the normal 
and tangential E— field components are comparable^ according to Eq. (J7]), i.e., \E X \ rs \E z \. 
This is in contrast to SP modes supported on metal-air interfaces that have /3' ~ ko, and 
have the normal electric field component as the dominant E— field component. For cases 
where the wavenumber of the SP supported on monolayer graphene satisfy (3' » ko, one 
should be able to make the approximation aj ~ (see Eq. (jHJ)), which remains accurate 
provided f3" « /3' and a" « a'. Substituting the approximation (3 ~ <x,- into Eq. (TTUj) 



yields a simple and physically intuitive solution 

(3 « iwe (ei + e 2 )/cr 
= J D(e 1 + e 2 )(u; 2 + ia;r- 1 )//i c , (12) 

where D = ireoh 2 /e 2 , and the substitution a ~ ^-mtra j^g i oeen m ade in the last step.— 
Separating Eq. f|T2|) into real and imaginary parts, we find 0' = D{e\ + e 2 )a; 2 //i c , and 
(3" = D(ei + e 2 )wr~ 1 //! c . This gives us some insight to the physical characteristics of the 
A mode. First, similar to the case of SPs supported on metal films, the A mode is better 
confined (larger (3') for higher values of ei or e 2 (as seen in Fig. 1(b)), accompanied with a 
proportional increase in absorption loss (larger /3"). Second, /?' varies as A~ 2 , decreasing at a 
rate much faster than j3" which varies as A -1 . Therefore, where the approximation is valid, n' 
varies as A -1 and n" is constant, in agreement with the dispersion curves of Fig. 1(b). Third, 
a higher doping level \i c decreases the confinement, which can be understood because the 
graphene becomes more like a perfect conductor with increased conductivity (see Eq. (12bl) ) 
and supports SP modes whose fields extends more into the dielectric media. This trend is 
shown in Fig. 2(a) with the case of isolated monolayer graphene (ei = e 2 = 1), where n! 
for different values of /i c is shown. It can be seen that the approximate solution Eq. (TT2|) 
is in excellent agreement with the exact dispersion relation fflOl) except near the plasmon 
resonance frequency. For antisymmetric A modes well- described by the approximation (1121) . 
the normalized propagation length L/X sp at which the field amplitude of the SP falls to 1/e 
of the initial value is^ 

L/\ sp « UT /2ix = cr/A , (13) 

which depends only on the operating wavelength and relaxation time r. As seen from 
the inset of Fig. 2(b), the approximate form cr/A remains highly accurate for operating 
wavelengths up to ~ 30 /im. The slight deviation for A > 30 /im is attributed to the 
inequality 0" << /3' in the approximation (Tl2|) . as the inequality is less accurate for longer 
wavelengths since as explained, f3' oc A~ 2 and f3" oc A -1 . The peak value of L/X sp is 
seen to increase with /i c according to the curve cr/A because a higher chemical potential 
has the effect of blue-shifting X p c , see Eq. ffTT]) . It is worthwhile noting that the ratio 
L/X sp = n' /2irn" can also act as a figure-of-merit (FOM) for ID plasmonic waveguides,— 
and is shown in Fig. 2(c) for 0.2 eV < fi c < 1.2 eV as a function of wavelength. However, 
in the case of graphene, L/X sp is virtually independent of important parameters like /i c , ei, 
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and e 2 (see Eq. ( lT3"j) ). As such we propose an alternative FOM, 

FOM = L/^\\7 P = Vn'/27rn" , (14) 

where the propagation length of the SP is normalized by the geometric mean of its effective 
wavelength and the free-space wavelength. The FOM as defined in Eq. (fi~4"|) for 0.2 eV < 
/i c < 1.2 eV as a function of wavelength is shown in Fig. 2(d). In contrast to the ratio L/X sp , 
the FOM is clearly distinct for different values of the chemical potential and wavelength. 
As such, we propose the newly defined FOM as a measure of the performance of graphene 
plasmon waveguides. For the graphene pair studied in the next section, this FOM will also 
be applied to quantify its performance as a SP waveguide. 

III. TUNING THE SP MODE WITH A PAIR OF PARALLEL GRAPHENE LAY- 
ERS 

Next, we demonstrate how one may tune the confinement and propagation characteristics 
of the SP mode by introducing an additional layer of graphene. Here, two monolayers of 
graphene (at x = 0, d) are separated by a thin dielectric medium (q) as shown in Fig. 3(a). 
Let us neglect the effects of absorption losses due to coupling with optical phonons (A op ~ 
6.3 /im) as described in Sect. I and consider the range of wavelengths A op < A < 50 /im. 
We are primarily concerned with graphene pairs whose spatial separation d is sufficiently 
large such that effects of interlayer hopping, which can be important for closely-spaced 
graphene pairs (few angstroms apart) or bilayer graphene, have no significant bearing on 
the SP dispersion.— For the following analysis, the smallest separation d we consider is 
eUin/A = 0.001, i.e. <i m i n ~ 6 nm. 

For a parallel pair of graphene monolayers separated by a distance d, the tangential 
H-field in each of the dielectric media can be expressed as 

NiH 1 exp(a\X + i(3z) x < , 

H -- < Ni(H a cosh adX + Hb smb. ad,x)exp(i(3z) < x < d , (15) 

NiH 2 exp(— a 2 (x — d) + i(3z) x > d , 

where N\ is a normalization constant similar to Nq in Eq. (jSJ). Matching the boundary 



conditions at x = and x — d, we obtain for the tangential if -field 



H a -(l + "^) Hl , 


(x = 0) 


(16a) 


iu 2 a 2 

H a cosh add + H b sinh a^<i = (1 H )-^2 , 

we 2 e 


(x = d) 


(16b) 


and for the tangential E z component, 








= 0) 


(17a) 


(if a sinh a rf ci + H b cosh a d c?) = H 2 , 


(x = d) 


(17b) 



where ad = (3 2 — q/cq. In Eq. (TT6|) . o"i and 02 are the conductivities of the graphene mono- 
layers at x = and x = d. The graphene conductivities o"i and a 2 , which are calculated with 
Eq. (CQ), may be tuned through their respective chemical potentials [i\ and fi 2 . Eliminating 
Hi,H 2 ,H a , and Hb from Eqs. (ITBT) and ( TT71) yields the dispersion relation 

tantieme? = - * 2 , (18) 

1 -f- 1 ii 2 

with 

r Od^ l^Oj 

<x,e d we^eo 

In the symmetric case where ei = e 2 and U\ = a 2 (/ii = /i2), the dispersion relation f fT8|) 
splits into two branches 

tanh— — = — — , (even, AS) (20a) 

2 I\ 

coth — = --j- , (odd, AA) (20b) 

2 Ti 

where Ti = T 2 due to symmetry. Let us take the case e\ = e 2 = e^ = 1, which is a config- 
uration that has also been studied in the context of Casimir forces^ and quasi-transverse 
electromagnetic waveguide modes^. Unless otherwise specified, it will be taken that the 
geometry is fully symmetric with \i\ = ii 2 = /i c , and e\ = e 2 = e^ = 1. 

For large separations d of the graphene pair, both Eqs. ffTBl) and (T201) reduce to the 
dispersion relation (fTOl) for the case of the I-G-I. For sufficiently small d such that the SP 
supported on individual graphene monolayer interacts with each other, the dispersion splits 
into an even and an odd mode. The typical field distribution of the if-field for these two 
modes are shown for d = 0.01A, A = 10 /im, and /i c = 0.4 eV as an example in Fig. 3(b). 
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The if —field displays a symmetric (antisymmetric) character across the gap for the even 
(odd) mode. Consistent with the analysis for monolayer graphene in Sect. II, the if —field is 
antisymmetric across each graphene layer. As such, let us label the even mode AS and the 
odd mode AA, where the first A in the labels refers to the antisymmetric ii-field distribution 
across each graphene monolayer. Additionally, let us note that the other two branches SA 
and SS (the first label S refers to the symmetric if- field across the graphene), can also be 
supported when the dielectric layer (e<j) is a thin metal film.— This case is verified with 
numerical simulations (not shown), where it was found that, similar to the case of the I-G-I, 
the dispersion with and without the graphene films differ only marginally. As such, we focus 
on the I-G-I-G-I structure. 

Figure 3(c) illustrates the splitting of the SP dispersion into the A A (thin blue curves) 
and AS (thick red curves) modes as the separation d is decreased, for the case A = 10 /im and 
/i c = 0.4 eV. Notice the difference in the vertical scales for n' and L/X sp . For sufficiently large 
d, both the AA and AS branches merge to the A mode supported by the monolayer graphene. 
As the separation d is decreased, the real part of the effective index n' increases (decreases) 
monotonically for the AS (AA) branch, corresponding to improved (degraded) confinement 
of the SP mode. On the other hand, the imaginary part of the effective index n" (inset 
of Fig. 3(c)) for both modes exhibits a nonmonotonic behavior. While this nonmonotonic 
behavior is also observed in the normalized propagation length L/X sp for the AA branch 
(dashed blue curve), it is absent for the AS branch (dashed red curve) however. This is 
because the SP wavelength decreases more quickly than the absorption loss increases as d is 
decreased. At very small separations d, or equivalently long wavelengths A, the AS branch 
would be strongly suppressed due to the exponential increase in absorption. 

For the AA mode, it is expected, within linear optics, that the amplitude of the if -field 
across the dielectric gap (e^) approaches asymptotically to a saw-tooth waveform with a 
near- zero amplitude as the separation d decreases (see Fig. 3(b)). This implies that the 
electromagnetic field in the gap is squeezed out of the gap into the two dielectric half-space 
as d decreases, explaining the decrease in the confinement (n') in Fig. 3(c). In the limit 
d — > 0, the AA mode tends to the A mode for the case of the bilayer graphene.— Here 
we estimate the conductivity of the bilayer graphene to be twice the conductivity of the 
monolayer, following the approximation Na for few layer graphene,— iM»M where N is the 
number of layers (N < 6). In Fig. 3(c), it is seen that n' for the AA branch approaches that 
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of the bilayer case, marked with a circle on the vertical axis. For \x c = 0.4 eV and A = 10 /im, 
the SP dispersion of the monolayer graphene satisfies the condition /3" « f3' (see Fig. 2(a)), 
and thus L/X sp for the bilayer case and monolayer case are almost equal (see inset Fig. 2(b)). 
These observations are all in agreement with the data for the AA mode represented by the 
blue curves in Fig. 3(c). For the AS mode, as long as the separation d is not too small 
(d > 0.01 A for this example), the increase in the confinement is accompanied with reduced 
absorption loss (n"). This concomitant improvement in confinement and absorption loss, 
which is contrary to the trade-off between confinement and loss associated with conventional 
SP modes,— is also observed for SPs supported in dielectric-loaded metallic waveguides^ 
that exhibit similar dispersion characteristics as the AS mode. 

The influence of the operating wavelength (A) and chemical potential (/x c ) on the dis- 
persion of the AA and AS modes are shown in Figs. 4 and 5, respectively. As observed, 
the behavior of the AS and AA branches share a number of similar characteristics with the 
case of the monolayer graphene investigated above. Let us first focus on the influence of 
the wavelength A (Fig. 4), where the chemical potential is taken to be \i c = 0.4 eV. For 
both branches, Figs. 4(a) and 4(b) show that n' is inversely proportional to A, whereas n" 
changes only weakly for different wavelengths (Figs. 4(c) and 4(d)). The result is an overall 
decrease in the normalized propagation length L/X sp as A increases. Additionally, because 
an increased wavelength effectively decreases the separation d, the two branches begin to 
merge for larger values of the normalized separation d/X. 

Figure 5 shows the influence of the chemical potential \x c on the dispersion of the two 
modes, where the wavelength A is taken to be 10 /im. Both n! (Figs. 5(a) and 5(b)) and 
n" (Figs. 5(c) and 5(d)) decreases as \x c increases, i.e., a higher chemical potential increases 
the conductivity of graphene, thereby reducing the confinement and the absorption loss. 
For a higher value of the chemical potential, because of the longer effective SP wavelength 
(X sp = X/n'), the separation between the graphene pair is effectively decreased, and thus the 
two branches start merging at larger separations d/X. A common feature of all the dispersion 
curves in Figs. 4 and 5 is that the change in n or n" becomes increasingly less profound as 
A or )jl c increases. These behaviors suggests an inversely proportional relationship between 
the wavenumber f3 of each SP mode (AA or AS) and A or /i c , in a similar fashion to Eq. (TT2|) 
for the monolayer case in Sect. II. From the above observations, the influence of e^ on 
the dispersion curves of the two modes may also be inferred. Since increasing q has the 
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same effect of reducing the effective SP wavelength X sp , the behavior for a greater dielectric 
constant of the separating medium (ea > 1) is qualitatively similar to a lower value of the 
chemical potential. As such, the AS and AA branches would merge for smaller separation d/X 
for a higher value of e^ (with all other parameters kept constant). Similar to a lower chemical 
potential, a higher e^ would result in a higher maximum value of n" . These behaviors for 
different values of e^ have been verified with calculations (not shown) . 

Next, let us compute the FOM as defined in Eq. Q14p for the AA and AS modes. We shall 
make a comparison between the FOM and the ratio L/X sp to see if the newly-defined FOM 
serves better to quantify the performance of the two modes for waveguiding operation. 
The two quantities L/X sp and FOM for the two modes as a function of the normalized 
separation d/X and wavelength A are shown in Fig. 6, with the chemical potential taken 
to be fi c = 0.8 eV. Except for a narrow region confined to d/X < 0.1 and A < 20 /xm, 
Fig. 6 (top panels) reveals that the ratio L/X sp varies only very weakly with d/X for a given 
wavelength. Furthermore, the ratio L/X sp fails to highlight the fact that the AS mode is 
extremely lossy for small separations d. By looking at the top panels of Fig. 6, one might 
also be misled to believe that the AA and AS modes are not useful for waveguiding operation 
for the long wavelengths in the terahertz regime. These deficiencies are circumvented with 
the corresponding FOM shown in the bottom panels of Fig. 6. First, let us note that unlike 
the ratio L/X sp , the FOM of the two modes for a given set of parameters lie in a similar 
range of values (see the scale on the colorbars). This is due to the factor \fn! in the FOM 
of Eq. Q14p instead of n! in the ratio L/X sp . Second, the FOM of the AA mode for very 
small separations d, or equivalently long wavelengths, can be relatively high. This is due 
primarily to the smaller absorption loss as the AA mode tends to the bilayer case. On the 
other hand, the AS mode provides a better FOM for intermediate values 0.05 < d/X < 0.2 
and 15 /mi < A < 45 /im. Thus for a properly chosen separation d, the parallel graphene 
pair offers at least one plasmon mode that is suitable for waveguide operation over the mid- 
infrared and terahertz spectrum. For a higher (lower) fi c , the absorption loss n" decreases 
(increases), confinement degrades (improves), and the general trend is an increase (decrease) 
in the FOM with the two branches starting to merge at larger (smaller) separations d/X. 

So far we have considered fully symmetric geometries. Yet another degree of freedom to 
tune the SP dispersion of the two branches is to control the chemical potential of one of 
the graphene monolayer while keeping the other one fixed. Evidently when \i\ ^ /12, the 
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symmetry is broken, and the fields of the modes no longer obey even or odd symmetry (see 
Figs. 7(a) and 7(b)). As such, let us define new branches AS' and AA' that stem from the 
AS and AA modes (respectively) for the symmetric case where /xi = /x 2 . As an example, let 
us take that /xi = 0.4 eV is fixed with /x 2 adjustable for an operating wavelength A = 30 /xm. 
For this case, the two modes for which /x x = /x 2 = 0.4 eV are the AA and AS modes, whose 
dispersion has been given in Fig. 4 (red dotted curves). The influence on the dispersion 
when /x 2 is varied so that \i\ and /x 2 are no longer equal is shown in Figs. 7(c) to 7(f). To aid 
the analysis, let us further categorized the dispersion behavior of the AS' and AA' branches 
for /x 2 < /xi and /x 2 > \i\. For the AS' (AA') branch, the dispersion is only weakly influenced 
for /x 2 > H\ (/x 2 < /xi), but changes visibly for /x 2 < fj,j (/x 2 > /xi). Figures 7(c) and 7(d) 
show that n' generally decreases for increasing /x 2 for both branches, consistent with earlier 
analysis on the AA and AS modes (see Figs. 4 and 5). As observed in Figs. 7(e) and 7(f), 
the dispersion of n" for the AA' and AS' branches exhibits a non-monotonic behavior, also 
similar to their symmetric counterparts. For the AS' branch with /x 2 < /xi and the AA' 
branch with /x 2 > /Xi, both n' and n" decrease for increasing /x 2 as expected. 

The dispersion behavior is however a little different for the weakly influenced case. In this 
case, for a qualitative explanation of the dispersion, let us take the AS' (AA') branch with 
A*2 > A*i (a*2 < Mi) to be a mode that has characteristics dominant of the AS (AA) mode but 
with a slight perturbation due to contribution from the AA (AS) mode, i.e., the two modes 
are mixed but with a dominant AS or A A character. For the AS' branch, n" is decreased 
very slightly for increasing /x 2 for d/X < 0.02. In the range 0.02 < d/X < 0.12, the AS' mode 
is modified with contributions from the more lossy AA mode, thus the absorption loss n" is 
increased even for /x 2 > /xi. For the AA' branch, n" is increased slightly for decreasing /x 2 for 
d/X < 0.04. In the range 0.04 < d/X < 0.12, the AA' mode is modified with contributions 
from the less lossy AS mode, thus decreasing the absorption loss n" for /x 2 < /Xi. 

We have shown the typical mode profile for the AS' and AA' modes for /x 2 ^ /xi in 
Figs. 7(a) and 7(b). In agreement with the analogy of the mixed mode, it is seen that the 
case of the weakly influenced dispersion (/x 2 = 0.8 eV for AS' given in dashed blue curve 
of Fig. 7(a) and , /x 2 = 0.3 eV for AA' given in solid red curve of Fig. 7(b)) displays a 
relatively stronger AS or AA behavior. On the other hand, the solid red curve in Fig. 7(a) 
(/x 2 = 0.2 eV) resembles that of the A mode (for monolayer graphene) with very little 
influence from the other graphene layer, and the dashed blue curve in Fig. 7(b) (/x 2 = 0.8 eV) 
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also resembles that of the A mode but slightly split within the graphene pair. A complete 
understanding of the physical mechanism leading to the weakly dispersive behavior of the 
AS' and AA' modes for the somewhat contrasting conditions described above is a future 
topic of investigation. The above analysis also shows that, while an asymmetric chemical 
potential generally increases the absorption loss for the AS' branch, it may serve to reduce 
the absorption loss for the AA' branch quite significantly. The corresponding FOM for 
the two branches are shown in Figs. 7(g) and 7(h), where it is seen that the dispersion 
characteristics described above are well-captured with the FOM. While the highest FOM is 
associated with the AA' branch for large values of /^ (/^2 > /-ti), the lowest FOM is observed 
for the AS' branch for small values of /i2 (/i2 < A*i)- 

IV. CONCLUSION 

In summary, for this article we have revisited the plasmon disperison in monolayer 
graphene, followed by a thorough investigation of the dispersion of the SP modes sup- 
ported in a parallel pair of graphene monolayers separated by a distance d. Due to the 
matching of boundary conditions across the interface of the monolayer graphene, it is found 
that the I-G-I combination supports only the A mode (antisymmetric magnetic field across 
the monolayer). To quantify the waveguiding performance of the SP mode, we found that a 
newly-defined FOM = L/ ^XX sp = y/n//2iin" is more optimal than the ratio L/\ sp because 
for a wide spectral range for which the approximation (fTSZj) is satisfied, the latter is almost 
independent of important parameters such as the chemical potential and dielectric constants 
of the neighboring media. 

In the case of the graphene parallel pair, the separation d is employed as an additional 
degree of freedom to tune the SP dispersion. For the fully symmetric geometry, the dispersion 
for the graphene pair splits into an even AS and odd AA branch for sufficiently small 
separations d. While the AS and AA modes may be characterized based on the symmetry 
of each mode, they share similar characteristics as the A mode supported in monolayer 
graphene. For instance, the influence of the wavelength and the chemical potential on the 
dispersion of the AS and AA modes are similar to that of the A mode in monolayer graphene, 
i.e. their wavenumbers are somewhat inversely proportional to the wavelength and chemical 
potential, similar to Eq (\V2\i . The non-monotonic behavior of the absorption loss (n") means 



15 



that there exist a certain range whereby the AS and AA modes are less absorptive compared 
to the monolayer case. By calculating the FOM, it is seen that the AA mode is optimal 
for small separations d or equivalently, long wavelengths A. On the other hand, the AS 
mode offers a better FOM than the AA mode for an intermediate range of separations and 
wavelengths where it is both less absorptive and better confined. For very small d however, 
the AS mode is strongly suppressed as its absorption increases exponentially. It is noted 
that the AA and AS branches here correspond to the optical and acoustic plasmon branches 
(respectively) in existing literature on theoretical investigation of the plasmon dispersion 
relation of the parallel graphene pair.— ^^ Varying the chemical potential of one of the 
graphene monolayers allows for yet another degree of freedom to tune the SP dispersion. 
Strictly speaking, the fields no longer possess symmetry in this case. However, we take the 
case /ii = Hi as reference, and consider branches AS' and A A' that stem from the AS and 
AA modes respectively. Our analysis shows that, under the right conditions, an asymmetric 
chemical potential can significantly reduce the absorption loss for the AA' branch. Our 
results demonstrate the potential to effectively tune the SP mode dispersion with a parallel 
pair of graphene monolayers to achieve well-confined and propagative SP modes in the mid- 
infrared and terahertz regime with graphene structures. This graphene plasmon waveguide 
platform is promising for mid-infrared and terahertz applications. 
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FIG. 1: SP dispersion characteristics for I-G-I structure, with £2 = 1. The geometry is illustrated in 
(a). The real part of the SP effective index n' for ei = 1 (green solid curve) and ei = 4 (red dotted 
curve) as a function of wavelength is shown in (b), inset: the corresponding imaginary part n" . In 
(b), the colored curves are results from Eq. (|10p . and black dashed curves depict calculations of the 
SP dispersion relation with the RPA2i shown here for comparison. Field components (E X ,E Z , H' = 
rjoH) with A taken to be 8 /xm are shown for (c) t\ = 1 and (d) e\ = 4. In (c) and (d), the normal 
electric field component E x is multiplied by a factor of 0.1 for the purpose of illustration. The 
chemical potential \i c is taken to be 0.8 eV for all calculations. 
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FIG. 2: SP dispersion characteristics for isolated monolayer graphene (ei = e% = 1). (a) n' 
and (b) FOM based on Eq. (j!4[) for several values of the chemical potential. The legend in (a) 
applies to (b). In (a), black dashed curves depict the approximate solution in Eq. (|12p . Inset of 
(b): normalized propagation length L/\ sp and the approximation cr/A (black dashed curve), (c) 
L/X sp and (d) FOM as a function of the chemical potential and wavelength. 
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FIG. 3: Splitting of the SP dispersion into the AA and AS modes with the geometrically symmetric 
I-G-I-G-I structure (ei = €2 = e<i = 1, /ii = ^2 = fJ-c), illustrated here for the case A = 10 /xm, 
and /i c = 0.4 eV. (a) Depicting the geometry, (b) Plot of the typical magnetic field component 
of the symmetric AS and antisymmetric AA mode with the separation d taken to be 0.01A. (c) 
Dispersion curves as a function of the normalized separation d/X for the AA (thick red curves) and 
AS (thin blue curves) mode. The circle on the left vertical axis in (c) indicates n' for the bilayer 
(I-2G-I) case. 
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FIG. 4: Influence of the operating wavelength on the dispersion of the AA and AS modes. Disper- 
sion curves for n' ((a) and (b)) and n" ((c) and (d)) of the AS and AA modes are shown in the 
top and bottom panels respectively, with the chemical potential fi c taken to be 0.4 eV. The legend 
in (c) applies to all other subfigures. 



21 



801 



(a) 



I'*— ■ 



v 40 k 



AS 



1.5 
1 



:(c) 



^ 



AS 



0.5 &-' 

oE 



40Y 



20 f, 





0.1 0.2 0.3 

dIX 



AA 



0.1 0.2 0.3 

dIX 



] r: 



(b) 

0.1 0.2 0.3 

d/X 



1.5 
1 

0.5 % 




A 




(d) 




AA 










*-^ 



0.1 0.2 0.3 

d/X 



H c = 0.2eV (i o = 0.3eV n c = 0.4 eV n c = 0.6eV |a c = 0.8 eV 

FIG. 5: Influence of the chemical potential on the dispersion of the AA and AS modes. Dispersion 
curves for n' ((a) and (b)) and n" ((c) and (d)) of the AS and AA modes are shown in the top 
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the bottom of the figure. 
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FIG. 6: Waveguiding performance of the AA and AS modes as a function of wavelength A and 
normalized separation d/X, with jj,\ = jj,2 taken to be 0.8 eV. Top panels: L/X sp , bottom panels: 
FOM. 
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FIG. 7: Influence on the SP dispersion when the chemical potential on one of the graphene sheets 
is varied (|Ui ^ 1x2). It is taken that [i\ = 0.4 eV is fixed, while [12 varies from 0.2 eV to 1.2 eV. 
Strictly speaking, the fields no longer possess symmetry for fi\ ^ fJ-2, as can be seen in (a) and 
(b). However, we take the case ^1 = ^2 = 0.4 eV as reference, and consider branches AS' and 
AA', which stem from the AS and AA modes respectively. Dispersion curves for n' ((c) and (d)), 
n" ((e) and (f)), and the corresponding FOM of Eq. (|14p ((g) and (h)) as a function of d/X for 
various values of \x-i for AS' and AA' are given in the top and bottom panels, respectively. The 
legend in (g), where the chemical potential \ii is specified in units of eV, applies to subfigures (c) 
- (h). The wavelength is taken to be A = 30 /mi. 
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